Mem. S.A.It. Vol. , 1 
©SAIt 2008 



„ono m/< ,^ 



Memorie deiia • 




Oegli S9" 



o 

(N 
00 



o 



- 1—1 

X 



Cosmological Numerical Simulations of Radio 
Relics in Galaxy Clusters: Insights for Future 

Observations 

Jack O. Burns 1 ' 2 and Samuel W. Skillman 1 ' 3 

1 Center for Astrophysics & Space Astronomy, University of Colorado at Boulder, Boulder, 
CO 80309 USA 

2 NASA Lunar Science Institute, Ames Research Center, Moffett Field, CA 94035 USA 

3 DOE Computational Science Graduate Fellow 
e-mail: jack.burns@cu.edu 



. Abstract. The acceleration of electrons at shock fronts is thought to be responsible for radio 

relics, extended radio features in the vicinity of merging galaxy clusters. By combining 
high resolution Adaptive Mesh Refinement Hydro/N-body cosmological simulations with 
, an accurate shock-finding algorithm and a model for electron acceleration, we calculate the 

^ |j expected synchrotron emission resulting from cosmological structure formation. From these 

simulations, we produce radio, SZE and X-ray images for a large sample of galaxy clusters 
along with radio luminosity functions and scaling relationships. We find that with upcoming 
radio arrays, we expect to see an abundance of radio emission associated with merger shocks 
in the intracluster medium. By producing observationally motivated statistics, we provide 
predictions that can be compared with observations to further our understanding of electron 
shock acceleration and kinematic structure of galaxy clusters. 

Key words. Cosmology: theory - Galaxies: clusters: intracluster medium - radiation mech- 
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1. Introduction electrons gyrating in ICM B-fields producing 

steep-spectrum, so-called "radio relics". 

Colliding galaxy clusters are insightful astro- 
physical plasma laboratories. Shocks produced Radio relics are rarely found in radio sur- 
during mergers heat the intracluster medium veys of galaxy clusters, possibly because of 
(ICM), assisting the gas to achieve roughly their steep radio spectrum and diffuse, low sur- 
hydrostatic equilibrium with the cluster grav- face brightness emission at cm wavelengths, 
itational potential well. Shocks also play a However, with new and improved radio ar- 
key role for the nonthermal component of the rays such as LOFAR, GMRT, and the EVLA, 
ICM. Shocks compress and amplify magnetic new high sensitivity observations may reveal 
fields. Shocks accelerate cosmic rays (CR) an abundance of cluster radio relics. To assist 
via a diffusive Fermi process. Thus, merger in searching for these relics, we have run high 
shocks are illuminated via the resulting syn- spatial dynamic range cosmological adaptive- 
chrotron radiation arising from relativistic CR mesh-refinement (AMR) simulations using the 
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Fig. 1. Projections of density (top), temperature (middle), and 1.4 GHz radio emission (bottom) for a rep- 
resentative sample of clusters from Enzo AMR simulations at z = 0. Each image is 4 lr l Mpc on a side. 
Peak resolution is 3.9 hT 1 kpc. 



Enzo code that produce ^2000 clusters in vari- 
ous st ages of merger evolution ( Ski llman et al. I 
120101) . Using new shock-identification tools 
and analyti cal models for diffusive shock ac- 
celeration dHoeft & Braggen1l2007l) . we have 
constructed radio maps of a large sample of 
simulated clusters to study the production and 
observational properties of relics in different 
merger states. D etails of these Enzo s imula- 
tions are found in Skillm an" et al. I (2010). 



2. Simulated Radio Relics 

Images of density, temperature, and radio 
emission for a small representative sample 
of clusters in a single projection along one 
axis are shown in Figure Q] Several inter- 



esting points are worth noting in comparing 
these images. First, there is a distinct differ- 
ence in morphologies between the density and 
radio images. The density (»X-ray emission) 
is center-filled whereas the radio is often edge- 
brightened. The curved radio arcs are illumi- 
nated bow shocks produced as two clusters 
pass between their cores. Such shocks are ef- 
fectively invisible on the density/X-ray maps 
and only partially visible on the temperature 
images. Thus, the radio relics light up im- 
portant evolutionary features in clusters (i.e., 
shocks) that are not apparent at X-ray energies. 

Second, in cases where the merger is 
largely along the plane of the sky, the radio 
relics fall on the edges of sharp temperature 
gradients. This is particularly apparent for the 
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Fig. 2. Profiles of X-ray flux, dimensionless 
Compton y SZE parameter (oc gas pressure), and ra- 
dio flux density for a merging cluster with a center- 
filled radio relic in the upper left-hand image of the 
lower panel in FigurefJJ 



cluster in the lower right of Figure [TJ This 
strong X-ray temperature/radio spatial correla- 
tion agrees with recent observations. 

Third, depending upon the projection, 
there is a wide variety of radio morpholo- 
gies. Although there is a preference for edge- 
brightened radio emission coinciding with 
shocks, there are a few clusters that demon- 
strate more diffuse, center-filled radio emis- 
sion (upper left-hand clusters in bottom panel 
of Figure [TJ. This corresponds to cluster pro- 
jections where the merger is largely along the 
line-of-sight. Such emission may qualitatively 
resemble "radio halos", but they do not demon- 
strate a strong correlation between X-ray and 
radio emission profiles seen for radio halos. In 
Figure [2] we show profiles of X-ray, Sunyaev- 
Zeldovich Effect (SZE), and radio emission 
for a cluster with apparently center-filled radio 
relic emission. For all the radio relics, the ra- 
dio profile shapes and slopes are inconsistent 
with those for the densities (and X-ray emis- 
sion). So, one should be able to distinguish be- 
tween real radio halos and projected relics via 
their X-ray/radio profiles, along with expected 
differences in spectral index. 



Fig. 3. Scaling relation between 1.4 GHz radio 
power and 0.2 - 12 keV X-ray luminosity for simu- 
lated clusters (solid line). Also shown (dashed line) 
is a line with the slope of the be st-fit observed scal- 
ing relation for radio relics from Ferettil d2002h . 

3. Scaling Relations 

In Fig. [3] we show the scaling relation be- 
tween radio power at 1 .4 GHz and soft X-ray 
luminosity for the most X-ray luminous sim- 
ulated clusters in our sample. Although there 
is a clear correlation between the radio and X- 
ray luminosities within r2oo ~ 0.8r,.„, a ;, there 
is considerable scatter in the relationship (fac- 
tor of « 10 4 in Piachz)- This scatter is real and 
represents different merger states for different 
clusters. High radio power clusters have suf- 
fered recent mergers whereas low power clus- 
ters have not experienced a merger in over a 
Gyr. To date, only the most radio luminous 
clusters have been observed and this possibly 
represents an observational bias that may be 
remedied with more sensitive, high bandwidth 
observations in the near future. 

There is also good agreement in the slopes 
of the scaling relation between observed and 
simulated clusters in Fig. [3] Although the ob- 
served sample from Feretti (2002) is small (9 
clusters), the general agreement between sim- 
ulations and observations is encouraging. 

4. Predicted Number of Radio Relics 

From our numerical simulations, we con- 
structed a radio luminosity function for radio 
relic clusters. That is, we calculated the cumu- 
lative number of clusters with P\achz greater 
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Fig. 4. Synthetic SZE/radio images. Contours are 
SZE flux and grey-scale is radio emission. 

than a given level, where the radio emission is 
the result of cluster merger shocks. Since radio 
power scales with cluste r mass in our simula - 
tions (Piachz K M i 2 ) (ISkillman et al. 11201 Oh . 
our computational volume produces clusters 
with M 20Q < 1O 15 M and P\agh z < 10 24 
W/Hz. To extend these limits to higher masses 
and radio powers, we extrapolated our clus- 
ter sample using the War ren et al. I (2006) mass 
function at < z < 1 . Within this cosmological 
volume of «26 (Gpc/h) 3 , our simulations and 
extrapolations predict that ^200 to 1000 ra- 
dio relic clusters with integrated (over an area 
with r = r 2 oo) P\agh z > 10 25 W/Hz will be 
present. With a factor of 10 improvement in 
the sensitivity of the EVLA at 1 .4 GHz, we ex- 
pect to see 10-100 times more galaxy clusters 
with radio relics. For an all-sky survey out to 
z ~ 0.5, there should be as 200 clusters with 

PlAGHz > 10 25 W/HZ. 

5. Shocks on SZE Images 

Fig.|4]shows an overlay of radio onto SZE im- 
ages for one cluster in our simulated sample of 
clusters. Note that the radio relic shocks cor- 
relate with sharp gradients on the SZE map. 
Also, there is extended SZE structure corre- 
sponding to the strongest relic emission. This 
may suggest that as observations continue to 



improve, SZE will complement radio images 
in mapping merger shocks. 

6. Conclusions 

Using a robust shock-finding algorithm and an 
analytical model for diffusive shock accelera- 
tion applied to a large-volume AMR cosmo- 
logical simulation, we have produced synthetic 
X-ray and radio images with characteristics 
similar to those observed in clusters with radio 
relics. The dual-arc radio morphologies and X- 
ray/radio scaling relations are good matches 
to observations. We predict that an increase in 
10-100 in the number of radio relic clusters 
detected with new and improved radio arrays 
(e.g., EVLA, LOFAR, and future lunar farside 
low frequency telescopes) should be possible. 
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